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m . Abstract 

We investigate the CP violating asymmetry, the forward backward asymmetry and the CP 
violating asymmetry in the forward-backward asymmetry for the inclusive B Xd £+^^decays 
■ for the £ = e, fi, t channels in the standard model. It is observed that these asymmetries are 

^ , quite sizable and B £+^^decays seem promising for investigating CP violation. 

o 

o ; 1 Introduction 

^ ; 

Qh, An efficient way in performing the precision test for the standard model (SM) is provided by the 

JlLiI flavor-changing neutral current (FCNC) processes since these are generated only through higher 

^ ' order loop effects in weak interaction. Among them, the inclusive B modes are 

prominent because of their relative cleanness compared to the pure hadronic decays. In the SM, 
B — > decays are dominated by the parton level processes b —>■ s{d)£^i^, which 

^ ■ occur through an intermediate u, c or t quarks. They can be described in term of an effective 

^ . Hamiltonian which contains the information about the short and long distance effects. 

The FCNC decays are also relevant to the CKM phenomenology; and b d£~^£~ modes 
are especially important in this respect. In case of the b —>■ s£^£~ decays, the matrix element 
receives a combination of various contributions from the intermediate t,coTu quarks with factors 
VtbVtl ~ A^, VcbV*, ~ and K^KTs ~ respectively, where A = sin (9c = 0.22. Since the 
last factor is extremely small compared to the other two we can neglect it and this reduces the 
unitarity relation for the CKM factors to the form VthV^* + VcbV*^ 0. Hence, the matrix element 
for the b — > s£^i^ decays involve only one independent CKM factor so that CP violation would 
not show up. On the other hand, as pointed out before fP.^, for b — > d£^l~ decay, all the CKM 
factors VtbV^^, VcbV*^ and VubV*^ are at the same order A'^ in the SM and the matrix element for 
these processes would have sizable interference terms, so as to induce a CP violating asymmetry 
between the decay rates of the reactions b — > d£'^£^ and b — > d£'^l~ . Therefore, b — > di^£^ 
decays seem to be suitable for establishing CP violation in B mesons. 

We note that the inclusive B — > Xsi~^£~ decays have been widely studied in the framework 
of the SM and its various extensions l3ll- fT9ll . As for B — > Xii£~^i~ modes, they were first 
considered within the SM in 1 1 1 and [T|. In ref. (2, together with the branching ratio, the CP 
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violating asymmetry for the B —^ X^^"*"^" decays has been studied including the long-distance 
(LD) effects, but only for £ = e mode. In Pl, a SM analysis for the forward-backward asymmetry 
is given again only for £ = e mode and neglecting the LD contributions. The general two Higgs 
doublet model contributions and minimal supersymmetric extension of the SM (MSSM) to the CP 
asymmetries were discussed in refs. [20] and [21 1, respectively. Ref. f2T| contains a comparative 
study of the CP asymmetries in the inclusive B X^i^ £^ and exclusive B — > ^ £~^£~ decays 
for £ = T only, by mainly focusing on the effects of the scalar interactions in the framework of 
the MSSM. Recently, CP violation in the polarized b — > d£~^i~ decay has been also investigated 
in the SM |22| and also in a general model independent way f23l. The aim of this work is to 
perform a quantitative analysis on the SM CP violation and the related observables, such as the 
forward-backward asymmetry and CP violation aysmmetry in the forward-backward asymmetry 
in the B —>■ X(ii~^£~ decays, some of which have already addressed in refs. l2l and l2Tl . 
as pointed out above. However, in this work we extend the investigation of the abovemensioned 
observables to consider all three lepton modes by mainly focusing on LD effects and also their 
dependence on the SM parameters p and r/. 

From the experimental side, the branching ratio {BR) of the B —>■ Xs£~^i~ decay has been 
also reported by the BELLE Collaboration |1T|, BR{B Xs£+i-) = ((6.1 ± 1.4)+{;^), which 
is very close to the value predicted by the SM |25 1, and may be used to put further constraint on 
the models beyond the SM. 

We organized the paper as follows: Following this brief introduction, in section |2l we first 
present the effective Hamiltonian. Then, we introduce the basic formulas of the double and dif- 
ferential decay rates, CP violation asymmetry, Acp, forward-backward asymmetry, ApB, and 
CP violating asymmetry in forward-backward asymmetry Acp{Apb) for B — > decay. 
Section|3]is devoted to the numerical analysis and discussion. 



2 The theoretical framework of B ^ X^i^i decays 

Inclusive decay rates of the heavy hadrons can be calculated in the heavy quark effective theory 
(HQET) 1 26 1 and the important result from this procedure is that the leading terms in l/rriq ex- 
pansion turn out to be the decay of a free quark, which can be calculated in the perturbative QCD; 
while the corrections to the partonic decay rate start with 1 / only. On the other hand, the pow- 
erful framework for both the inclusive and the exclusive modes into which the perturbative QCD 
corrections to the physical decay amplitude are incorporated in a systematic way is the effective 
Hamiltonian method. In this approach, heavy degrees of freedom, namely t quark and bosons 
in the present case, are integrated out. The procedure is to take into account the QCD corrections 
through matching the full theory with the effective low energy one at the high scale /x = mw and 
evaluating the Wilson coefficients from mw down to the lower scale ij. ~ 0{mb). The effective 
Hamiltonian obtained in this way for the process b —>■ d£^£~, is given by lfT4l . I27l - ll30l : 

We// = ^ Vm!^ X] C,{p) O^ifi) - A„{Ci(/i)[On/") - Oi(/i)] 

+C2{p)[OUf^) - 02{p)]}] (1) 



where 



= -T^, (2) 
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using the unitarity of the CKM matrix i.e. VtbV/^ + VubV*^^ = —VcbV*^. The explicit forms of 
the operators Oi can be found in refs. l27l l28l . In Eq.Q, Cj(^) are the Wilson coefficients 
calculated at a renormalization point ^ and their evolution from the higher scale /x = down 
to the low-energy scale /i = nib is described by the renormalization group equation. For Cj^^{fi) 
this calculation is performed in refs. llBTIlB? ! in next to leading order. The value of Cio{mb) to the 
leading logarithmic approximation can be found e.g. in 1,27. .30,1 . We here present the expression 
for Cg{fi) which contains the terms responsible for the CP violation in i? — > Xdi~^£~ decay. It 
has a perturbative part and a part coming from long distance (LD) effects due to conversion of the 
real cc into lepton pair £~^£~ : 

Cg"(.f^) = Cr*(M) + Yresonis) , (3) 

where 

+ Xu{3Ci + C2)] - s) (4C3(/i) + 4C4(/u) + 3C5(/i) + Ce{fi)) 

- h{0,s) [Csifi) + 3C4(/i) + A„(6Ci(/i) + 2C2(m))] (4) 

+ l{3C3{fi) + Ciifi)+3C5{fi) + Ceifi)) , 

and 

X [(3C7i(;u) + C2{ii) + 3C3(;u) + C4(/x) + 3C5(/i) + C^{ix)) 

+ A„(3Ci(/x) + C2(a^))]. (5) 

In Eq.©, s = jrr?^ where q is the momentum transfer, u = ^ and the functions s) arise 
from one loop contributions of the four-quark operators 0\ — 0% and are given by 

, / N 8 , mf, 8 , 8 4 

n(u,s) = — m muH 1 — y (6) 

^ ' ' 9 /i 9 27 9 



2 -, ,„ In 

--(2 + y)|l-y|i/2 



ivr ) , for y = ^ < 1 



2 



s 



2 arctan for y = — > 1, 



, / X 8 8 4 , 4 

M0,5 = ---ln^--lns + -i7r. (7) 
27 9 /i 9 9 

The phenomenological parameter k in Eq. (|5ll is taken as 2.3 (see e.g. L33J ). 

The next step is to calculate the matrix element of the B Xdt^l~ decay. Neglecting the 
mass of the d quark, the effective short distance Hamiltonian in Eq.Q leads to the following QCD 
corrected matrix element: 

M = ^I^VM\cl^f{mb) - 75)6^V^ + C7io(mb) d-f^{l - 75)6^^75^ 

- 2C7^^^(m6)^(ila^.(z"(l+75)&^V^|. (8) 
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Since the initial and final state polarizations are not measured, we must average over the initial 
spins and sum over the final ones, that leads to the following double differential decay rate 



ds dz 



a 



* |2 
td\ 



+ 12 (l + -) Re{C^^^Cg^^*) + 6f Re(C7ioC. 
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3 

+ 2 
+ 6 
3 

+ 2 



a 



1 + s) - (1 - s)v'^z'^ + 4t 

1 ]v^z^ H 

s/ s 



eff\2 



1 

1 + - 

s 



2^2 



(1 +s) - (1 - s)v^z 



At 



\Cio\'} 



(9) 



where v = ^/l — At/s, t = mj/ml and z = cos 6, where 9 is the angle between the momentum 
of the B-meson and that of in the center of mass frame of the dileptons i^i'^. In Eq. Q, 



Glml 
1927r3 



|^c6|'/(n)A;(n), 



where 



k{u) 



1 - 8n + 8n*^ - ti^ - 2AuHn{u) 



1 



2a^(mb) 



vr 



31 

T 



)(i-^c) + 



(10) 

(11) 

(12) 



are the phase space factor and the QCD corrections to the semi-leptonic decay rate, respectively, 
which is used to normalize the decay rate of i? — > Xdi'^l~ to remove the uncertainties in the 
value of TTife. 

After integrating the double differential decay rate in Eq.® over the angle variable, we find 



(1-.) 



* |2 



4t . , X 
1--A s , 
s 



(13) 



where 

Ms) 



(s + 2s2 + 2t- 8st) 



Cio\^ + -{2 + s){s + 2t)\Cff\^ 
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+ —{s + 2t)Re{C^^^Cg 



(2 + «)(l + f)|C9^^^P 

(14) 



We start with calculating the CP asymmetry Aqp between the B —>■ Xdi~^£ and the conju- 
gated one B Xcii~^£~ , which is defined as 



where 



Acp{s) 

dr dr{B xj+i-) 



dr _ dr 

ds ds 
ds ds 



ds 



ds 



dr _ dr{B^Xdi+t 

ds 



ds 



(15) 



(16) 



Since in the SM only Cg^^ contains imaginary part, representing Cg^^ symbolically as 



a 



eff 



(17) 
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and further substituting A — > A* for the conjugated process B —>■ X^l'^l , one can easily obtain 
lil 

, , , -2Im(A„)S 

where 

E = (l + -)[(l + 2s)Im(Ci*6)+6C7^^^Im(e2)]Im(AJ. (19) 

For completeness, we next consider the forward-backward asymmetry, ApB, vaB ^ 
which is another physical quantity that may be useful to test the theoretical models. Using the def- 
inition of differential ApBis) 



Jo "-^dsdz J-1 
JO "'^dsdz ^ J-1 "-^dsdz 



A / „\ JO ""^ dsdz J — l"'"dsdz /'OA\ 

Afb[s) = Pi , ^2r , pO , d^T ' ^^^> 



we find 



Afb{s) = A^Re[Cio(2C7^^^ + sC9^^^*)], (21) 

which agrees with the result given by ref. 13, but not by l2Tl . 

We have also a CP violating asymmetry m ApB, Acp{Apb), in i? — > £+£~decay. Since 
in the limit of CP conservation, one expects App = —App f2"34l, where App and ^i?^ are the 
forward-backward asymmetries in the particle and antiparticle channels, respectively, Acp{Apb) 
is defined as 

Acp{Apb) = App + App . (22) 
Here, App can be obtained by the replacement, 

Ct^f{\u)^Cl"{\u^\l). (23) 

Using Eqs. d^ we can find 

4 ,A N 6 V Im(A„) 

= A(A + 4Im(A„) 

• 2S {2C^^^^ + s(Re(ei) + Re(e2) Re(A,) - Im(6) Im(A„))) - s A Im(e2)] , 

(24) 

which is sUghtly different from the one given by ref. 1^ . 



3 Numerical analysis and discussion 

In this section, we present results of our calculations related to i? ^ decays, for two 

different sets of the Wolfenstein parameters. For this we first give the Wolfenstein parametrization 
(351 of the CKM factor in Eq.© 



p{ l - p)-rf 
(l-p)2+r?2 



A. = , ' + O(A^) , (25) 
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and also 



A2[(1-p)2 + ^2j + 0(A4). (26) 



The updated fitted values for the parameters p and r] are given in ref. l36l as 

p = 0.22 ±0.07 (0.25 ±0.07) 

f] = 0.34 ± 0.04 (0.34 ± 0.04) (27) 

with (without) including the chiral logarithms uncertainties. In our numerical analysis, we have 
used (/), iij) = (0.15; 0.30) and (0.32; 0.38), which are the lower and higher allowed values 
of the parameters given in Eq. (l27l above, and present the dependence of the Acp, ApB and 
Acp{Afb) on the dimensionles photon energy s for the B X(ii^£^{i = e, /i, r) decays in 
Figs. (fTl6l . 

We have also evaluated the average values of CP asymmetry < Acp >, forward-backward 
asymmetry < App > and CP asymmetry in the forward-backward asymmetry < Acp{Apb) > 
in i? — > decay for the above sets of pai^ameters (/?, r]), and our results are displayed in 

Table 1 and 2 without and with including the long distance effects, respectively. 

The input parameters and the initial values of the Wilson coefficients we used in our numerical 
analysis are as follows: 

mB = 5.28 GeV , = 4.8 GeV , = 1.4 GeV , mt = 175 GeV , 

me = 0.511 MeF , = 1.777 GeV, = 0.105 GeF, , 

BR{B X^eUe) = 10.4% , a = 1/129 , mw = 80.4 GeF , mz = 91.1 GeF 

Gi = -0.245, G2 = 1.107, G3 = 0.011, G4 = -0.026, G5 = 0.007, 

Ge = -0.0314, C^^^ = -0.315, Gg = 4.220, Gio = -4.619. (28) 

In our numerical analysis, we take into account five possible resonances for the LD effects 
coming from the reaction b —>■ dtpi ^ di~^£~ , where i = 1, 5 and divide the integration region 
into two parts for £ = t: {2m^/mB)'^ < s < {{m^^ — 0.02)/mB)^ and ((m^^ ± 0.02)/mB)^ < 
s < 1, where m^^ = 3.097 GeV is the mass of the first resonance. As for £ = e and /x modes, 
the integration region is divided into three parts : {2m^/mB)'^ < s < {{m^^ — 0.02)/mB)^, 
((m^j ± 0.02)/mB)2 < s < ((m^^ - 0.02)/mij)2 and ((m^^ + Qm)/mBf < s < 1, where 
m^2 = 3.686 GeV is the mass of the second resonance. 

For reference, we present our SM predictions with long distance effects 

BR{B^Xd£^£') = (3.01,2.61,0.11) X 10"^ (29) 

for i = e, p, T, respectively, with {p; rj) = (0.30; 0.34), which is in agreement with the results of 
ref.ni. 

In Fig.Q and Fig.©, we present the dependence of Acp on the dimensionless photon energy 
s, for B Xci£~^£~deca.y for the Wolfenstein parameters (p; r/) = (0.15; 0.30) and (p; ry) = 
(0.32; 0.38), respectively. The three distinct lepton modes £ = e, p, t we represented by the 
dashed, dotted and solid curves, respectively. We observe that the Acp for £ = e, p cases almost 
coincide, reaching up to 25 % for the larger values of s. The Acp for ^ = r mode exceeds the 
values of the other modes and reaches 40 %. We also observe from Tables 1 and 2 that including 
the LD effects in calculating < Acp > does not change the results for £ = e, p modes, while 
£ = T mode, it is quite sizable, 8 — 36%, depending on the sets of the parameters used for (p; t]). 

The s dependence of App for the B —^ X^ £^£'' {£ = e, p, r) decays are plotted in Figs.(|3ll 
and © for {p; rj) = (0.15; 0.30) and (p; r/) = (0.32; 0.38), respectively. We see that App is 
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(p; v) 


< Acp > 
I = e £ = At e = T 


< ApB > 

I = e e = fi e = T 


< AcpiApB » 

e = e e = e = T 


(0.15; 0.30) 


0.030 0.036 0.134 


-0.124 -0.151 -0.182 


-0.009 -0.009 0.001 


(0.32; 0.38) 


0.051 0.061 0.169 


-0.129 -0.156 -0.180 


-0.015 -0.015 0.002 



Table 1: The average values of Acp, Aps and Acp{Apb) in i? — > X^i'^i for the three distinct 
lepton modes without including the long distance effects. 





< Acp > 

e = e e = fj. i = T 


< ApB > 

e = e e = fi e = T 


< AcpiApB >) 

I = e e = fj. e = T 


(0.15; 0.30) 


0.032 0.036 0.144 


-0.119 -0.139 -0.157 


-0.017 -0.017 -0.004 


(0.32; 0.38) 


0.051 0.059 0.230 


-0.125 -0.140 -0.150 


-0.031 -0.030 -0.009 



Table 2: The same as Table (dJ, but including the long distance effects. 



negative for almost all values of s, except in the resonance and very small-s regions. < App > 
takes the values between —(12 — 15)% depending on the sets of the parameters used for {p; t]) . 
The LD effects on < App > are about 10%, but in reverse manner, decreasing its magnitude in 
comparison to the values without LD contributions. 

We present the dependence of the AcpiApp) of B ^ Xdi^£~decay on s in Fig.© and 
Fig-®, again for two different sets of the Wolfenstein parameters. As for Acp, Acp{Apb) for 
£ = e, and £ = modes almost coincide. We see that Ac p (App) is all negative except in a 
very small region for the intermediate values of s for £ = e, fj. cases. LD effects seem to be quite 
significant for < Acp{App) >, enhancing its value twice (four times) for £ = e, p {£ = r) 
modes. To see this LD contributions more closely, we present the < Ac p (App) > for different 
regions of s in Table (|3ll and (|4]i, for (p; r/) = (0.15; 0.30) and (/); r/) = (0.32; 0.38), respectively. 
We see that for the light lepton modes, £ = e, fi, AcpiApp) is more sizable in the high-dilepton 
mass region of s, {{m^^ + 0.02)/mpy < s < 1. However, for £ = t, the contribution from 
the high-dilepton mass region of s is negligible and the contribution to < Ac p (App) > comes 
effectively from the low-dilepton mass region, {2mi/mpY < s < (("^Vi ~ 0.02)/mB)^ and 
amounts to —1%. 

As a conclusion we can say that there is a significant Acp and Acp{Apb) for the B — > 
decay, although the branching ratios predicted for these channels are relatively small 
because of CKM suppression. So, B X^^+^^decays seem promising for investigating CP 



£ 


SD 
contribution 


{2mi/mpY < s < 
{{m^, - 0.02)/mB)2 


{{m^, + 0.02)/mB)2 < s 
< ((m^, - 0.02)/mpy 


((m^, + 0.02)/mB)2 
< s < 1 


SD+LD 
contribution 


e 


-0.92 


-0.29 


-0.25 


-1.20 


-1.78 


/i 


-0.91 


-0.29 


-0.25 


-1.20 


-1.78 


r 


-0.11 


-0.42 


3.10 X 10-3 




-0.42 



Table 3: The SM predictions for the average CP- violating asymmetry in the forward-backward asym- 
metry < Acp{Apb) > xlO^^ for different regions of the dimensionless photon energy s with 
{p;r]) = {0A5; 0.30). 



7 



I 


SD 
contribution 


{2mi/mBY < s < 
((m^, - 0.02)/mB)2 


((m^, + 0.02)/mB)^ < s 
< ((m^, - 0.02)/mB)2 


((mw,2 + 0.02)/mB)2 
<s<\ 


SD+LD 
contribution 


e 


-1.59 


-0.51 


-0.43 


-2.15 


-3.10 




-1.57 


-0.51 


-0.43 


-2.15 


-3.09 


r 


0.20 


-0.94 


3.30 X 10"=^ 




-0.94 



Table 4: Same as Table ©, but with (p; r]) = (0.32; 0.38). 



violation. 
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Figure 1: Acp for B X^t 
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Figure 2: The same as Fig.© but for the Wolfenstein parameters (p, r/) = (0.32; 0.3^ 
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Figure 3: Aps for B decay for the Wolfenstein parameters (p, rj) = (0.15; 0.30). The 

three distinct lepton modes i = e, /x, r are represented by the dashed, dotted and solid curves, 
respectively. 
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Figure 4: The same as Fig.© but for the Wolfenstein parameters (p, rj) = (0.32; 0.38) 
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Figure 5: Acp{Afb) for S decay for the Wolfenstein parameters (p, rj) = (0.15; 0.30) 

The three distinct lepton modes i = e, p, r are represented by the dashed, dotted and solid curves 
respectively. 
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Figure 6: The same as Fig.® but for the Wolfenstein parameters (p, t]) = (0.32; 0.3^ 
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